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ABSTRACT 


In  this  paper  boundary  value  problems  on  infinite  intervals  are 
treated.  There  is  a  theory  for  problems  of  this  kind  which  requires  the 
fundamental  matrix  of  the  system  of  differential  equations  to  have  certain 
decay  properties  near  infinity.  The  aim  of  this  paper  is  to  establish  a 
theory  which  holds  under  weaker  and  more  realistic  assumptions.  The  analysis 
for  linear  problems  is  done  by  determining  the  fundamental  matrix  of  the 
system  of  differential  equations  asymptotically.  For  inhomogeneous  problems  a 
suitable  particular  solution  having  a  'nice*  asymptotic  behaviour  is  chosen 
and  so  global  existence  and  uniqueness  theorems  are  established  in  the  linear 
case.  The  asymptotic  behaviour  of  this  solution  follows  immediately.  Non¬ 
linear  problems  are  treated  by  using  perturbation  techniques  meaning 
linearization  near  infinity  and  by  using  the  methods  for  the  linear  case  . 
Moreover,  some  practical  problems  from  fluid  dynamics  and  thermodynamics  are 
dealt  with  and  they  illustrate  the  power  of  the  asymptotic  methods  used. 
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SIGNIFICANCE  AND  EXPLANATION 


Boundary  value  problems  on  infinite  intervals  occur  in  many  areas  of 
physics,  for  example  in  fluid  mechanics  when  similarity  solutions  of  the 
Navier-Stokes  equations  describing  the  flow  over  an  infinite  medium  are 
sought.  Problems  which  are  dealt  with  in  this  paper  have  the  following 
form.  We  look  for  a  solution  of  a  system  of  differential  equations  over  an 
infinite  interval  which  is  continuous  at  infinity  and  which  fulfills  certain 
boundary  conditions  at  a  finite  point  and  at  infinity.  The  following 
questions  arise  immediately:  which  condition  on  the  differential  equation  and 
on  the  boundary  conditions  assure  the  existence  and  the  uniqueness  of  a 
solution  and  how  does  this  solution  depend  on  the  data .  There  is  a  well  known 
theory  for  problems  of  this  kind,  but  it  can  only  be  applied  to  problems  where 
the  fundamental  matrix  of  the  system  of  differential  equations  has  certain 
convergence  properties  near  infinity.  However  this  assumption  is  not 
fulfilled  for  many  practical  problems.  This  paper  answers  the  above  named 
question  under  very  weak  assumptions  on  the  problems.  Moreover  some  fluid 
dynamical  problems  illustrating  the  power  of  the  theory  developed  are  dealt 
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ANALYSIS  OF  BOUNDARY  VALUE  PROBLEMS 
ON  INFINITE  INTERVALS 


PETER  A.  HARROW ICH 

1 .  Introduction. 

This  paper  is  concerned  with  the  analysis  of  boundary  value  problems  on  infinite 
intervals  posed  in  the  following  way 

(1.1)  y'  =  t0f ( t,y) ,  1  <  t  <  ”  ,  a  c  NQ 

(1.2)  y  e  C([1,<»I):  <”*=>  y  £  C({1,“))  and  lim  y(t)  -  y { °» )  exists 

(1.3)  b(y(  1 )  ,y<°>) )  =  0 
where  y  is  an  n-vector,  f  and  b  are  nonlinear  mappings. 

Equation  (1.1)  has  a  singularity  of  the  second  kind  at  t  =  “  of  rank  a  +  1 , 
because  we  assume  that  f  is  continuous  in  (•,y(“)).  The  goal  is  to  establish  existence 
and  uniqueness  theorems  for  very  general  f's  and  b's  ,  to  gain  information  on  the 
behaviour  of  y  for  large  t  and  -  in  the  linear  inhomogenous  case  -  to  investigate  the 
dependence  of  y  on  the  boundary  data  and  the  inhomogenity. 

Problems  of  this  kind  frequently  occur  in  fluid  dynamics  when  stationary  similarity 
solutions  of  the  Navier-Stokes  equations  for  certain  flow-constellations  are  sought  (see 
for  example  McLeod  (1969);  Markowich  (1980);  Lentini  and  Keller  (1980),  Cohen,  Fokas  and 
Lager strom  (1978)). 

For  application  in  other  areas  of  physics  see  for  example  Lentini  (1978). 

Much  analytical  work  has  been  done  on  singular  boundary  value  problems  of  the  second 

d  f 

kind.  F.  de  Hoog  and  Richard  Weiss  ( 1980a  ,b)  investigated  the  case  where  («,,y(“‘)l  has 
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no  eigenvalue  on  the  imaginary  axis  by  linearizing  f  around  y(“)  and  evaluating  at 
t  “  so  getting  the  constant  coefficient  problem 

(1.5)  z'  -  tafy(“,y(“) )z 

and  by  employing  perturbation  techniques  which  are  based  on  estimates  derived  for  a  certain 

particular  solution  of  linear  inhomogenous  problems.  They  could  establish  uniqueness  and 

existence  theorems.  M.  Lentini  and  H.  B.  Keller  (1980)  extended  this  approach  neglecting 

the  assumption  on  the  eigenvalues  but  they  required  that  the  projection  of  f  (t,y(">) 

onto  the  direct  sum  of  invariant  subspaces  of  f^(“,y(“))  which  correspond  to  an  imaginary 

-  ( a+ 1 ) t—  e 

eigenvalue,  converge,  at  least  as  t  ,  where  r  is  the  largest  dimension  of  these 

subspaces  and  e  >  0  .  It  turns  out  that  this  assumption  is  crucial  and  the 
perturbation  approach  breaks  down  if  f^fm.yt*))  has  eigenvalues  with  real  part  zero  and 
if  the  convergence  requirement  is  neglected.  However  many  practical  problems  do  not 
fulfill  this  convergence-requirement  in  the  presence  of  Imaginary  eigenvalues  (see  for 
example  Cohen,  Fokas  and  Lagerstrom  (1978)  and  M.  Lentini  and  V.  Pereyra  (1977)  and 
therefore  a  more  general  approach  is  necessary. 

In  this  paper  asymptotic  series  are  used  in  order  to  determine  asymptotically 
fundamental  solution  matrices  of  linear  systems  of  the  form: 

(1.6)  y*  “  t°A(t)y  t  >  5  . 

The  basic  assumption  is  that  A  is  analytical  in  for  some  6  >  1  ,  so  that 

°°  e  i 

(1.7)  A(t)  »  f  A,  t_i  where  A  =  7-  lim  A(J)  . 

i  i  il  .  .  i  x 

i*0  X+0+  dx 

Then  a  fundamental  matrix  ♦ ( t )  of  (1.6)  can  be  calculated  as  an  asmptotic  (formal)  log- 
exponential  series  from  the  coefficients  A^  by  a  recursive  algorithm  (see  Coddington  and 
Levinson  (1955)  and  Wasow  (1965)). 

Assumption  (1.7)  can  be  weakened  so  that  only  a  finite  but  large  enough  number  of 
these  derivatives  exist.  Equation  (1.1)  is  treated  by  linearization  around  y(“) 
obtaining  the  variable  coefficient  problem 

(1.8)  z'  =  tafy(t,y(»))z 

and  again  by  employing  similar  perturbation  techniques.  The  advantages  o'  the  formal- 
series  approach  is  twofold.  Firstly  no  restrictions  (except  (1.7))  have  to  be  made 
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concerning  the  convergence  behaviour  of  fy(t,y{») ),  secondly  the  asymptotic  behaviour  of 
the  (basic)  solutions  is  obtained  directly.  The  asymptotic  behaviour  is  crucial  for  the 
determination  of  appropriate  numerical  procedures  for  problsms  (1.1),  (1.2),  (1.3).  (see 
Lentini  and  Keller  (I960),  Markovich  (1980)  and  de  Hoog  and  Weiss  (1980b)). 

This  paper  is  organised  in  the  following  way,  in  chapter  1  some  remarks  are  made  on 
linear  lnhomogenous  constant  coefficient  problems  (see  Lentini  and  Keller  (1980)),  in 
Chapter  2  the  case  where  A(")  has  distinct  eigenvalues  is  treated,  in  Chapter  4  we  admit  a 
general  Jordan  canonical  form  of  A(“)  and  ii  hapter  5  we  get  to  nonlinear  problems  of 
the  form  (1.1),  (1.2),  (1.3).  Chapter  6  is  concerned  with  practical  examples  which 
illustrate  the  power  of  the  used  asymptotic  methods. 


2.  Linger  probing  with  constant  coefficient*. 

We  consider  problems  of  the  form 

(2.1)  y'  -  t°Ay  «  taf(t)  1  <  t  <  -  ,  a  e  R  ,  a  >  -1 

(2.2)  y  e  C<n,“J) 

(2.3)  B^O)  +  B.prf-)  -  8 
where  the  n  x  n-Matrix  A  *  0  . 

At  first  we  transform  A  to  its  Jordan  canonical  form  J 

(2.4)  A  -  E  J  E-1 
and  substitute 

(2.4a)  u  -  E_1y  . 

So  we  get  the  new  problem 

(2.5)  u'  -  ta  J(t)u  =  taE_1f(t) 

(2.6)  u  e  C( (1,-1 ) 

(2.7)  BE  u(1)  +  BE  u{»)  =  6  . 

1  00 

Without  loss  of  generality  we  can  assume  that  J  has  the  block  diagonal  form 

+  0 

(2.8)  j  -  diag(J  ,  J  ,  J  ) 

where  the  real  part  of  the  eigenvalues  of  J+  are  positive.  The  real  parts  or  the 
eigenvalues  of  J°  are  equal  to  zero  and  the  real  parts  of  the  eigenvalues  of  J-  are 
negative.  This  structure  can  always  be  obtained  by  reordering  the  columns  of  E  .  Let  the 
dimensions  of  these  three  matrices  be  r+,  rg  resp.  r_. 

The  diagonal  projections  D+,  Dg,  D_  are  obtained  by  taking  the  main-diagonal  of  J 
and  by  replacing  every  eigenvalue  with  positive,  zero  resp.  negative  real  part  by  1  and 
all  others  by  zero  so  that 

(2.9)  I  -  D+  *  Dg  +  D_ 
holds. 

Furthermore  let  be  the  projection  onto  the  direct  sum  of  elgenspaces  of  J 

associated  with  zero  eigenvalues,  which  is  obtained  by  replacing  by  zero  every  (diaqonal) 
element  of  Dg  which  is  not  associated  with  the  first  column  of  a  Jordan  block  of  J 
belonging  to  a  z»ro  eigenvalue. 
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Let  the  number  of  nonzero  columns  of  DQ  which  equels  the  geometrical  multiplicity  of 
the  eigenvalue  zero,  be  rQ  .  The  general  solution  of  the  homogenous  problem  (2.5),  (2.6) 
iss 

a+1 

(2.10)  u^t)  “  +  “  «xp(j  ~j"  )  <5q  +  D_)5  ,  {  t  c”. 

In  order  to  solve  the  inhomogenous  problem  (2.5)  we  look  for  a  particular  solution 
up  e  C( [1,»J ). 

F.  de  Hoog  and  R.  Weiss  ( 1980a, b)  and  M.  Lentini  and  H.  B.  Keller  (1980)  suggested  the 
following  choice: 

u  (t)  -  ( Hf ) ( t)  -  $<t)  !t  D+*-1(s)E'1f(s)s  ds  + 

P  ao 

(2#11)  +  ♦< t)  /fc  Dq*  1 ( s )E  1 f ( s ) s^ds  + 

+  $(t)  ft  D  $  ^ ( s )E  1f(s)sads 
6 


with  6  c  [1,“)  . 

We  denote  the  three  terms  on  the  right  hand  side  of  (2.11)  by 

u  “  H  f,  u  *  H  f  resp.  u  “  H  f.  F.  de  Hoog  and  Richard  Weiss  (1980a)  showed 
P+  +  P0  °  P.  " 

that  u  and  u  are  in  C((1»“!)  if  D  E_1f  and  D.E-1f  are  in  C {  [  1  ,“* ]  )  and  that 
P+  P_  " 

J (D  +  D  )u  (“>  -  -(D  +  D  )E_1f(“)  holds. 

+  -  p  + 

M.  Lentini  and  H.  B.  Keller  (1980)  showed  that 

lu  ( t )  i  -  0(t"E)  if  ID  E_1fl  =■  o(t'(a+1  ,r“E) 

P0 

where  e  >  0  and  r  is  the  maximal  dimension  of  the  invariant  subspaces  of  J  associated 
with  eigenvalues  on  the  imaginary  axis.  Therefore  the  operator  H  operates  on  the  space 
of  all  functions  f  ,  which  fulfill: 

(2.12)  feC([1,-]>  and  Dqe‘ 1 f ( t)  =  FQ ( t ) t* <a+1 ’ r_E 

with  F.  c  C.  where  C.([1,—  )>  is  the  space  of  functions  which  are  continuous  on 

Ob  b 

[ 1  /•)  and  bounded  as  t  ♦  •  • 
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Inserting  the  general  solution  of  (2.5),  (2.6)  into  (2.7)  we  get 


[(B  E  +  BJSlD  +  B  Eexp(-^r)D_]C  - 

(2.13) 

-  0  -  [B,EU  (1)  B  Eu  (»)]. 

Ip  “  P 

Therefore  we  conclude 

Thereom  2.1.  The  problem  (2.1),  (2.2),  (2.3)  has  a  unique  solution  y  for 


*  rQ  +  r- 

sll  f  which  fulfill  (2.12)  and  Be  R  ,  if  and  only  if 


(2.14)  rankt^E  +  BJSJDg  +  B^xp^JoJ  -  tQ  +  r_ 
where  B,  and  B  are  (r„  +  r  )  x  n  matrices. 

|  OO  0  — 

In  this  case  y  depends  continuously  (in  the  norm  I y 1  :» 

1 1  ,WJ 

* 

max  ly(t)l)  on  the  data  B  ,  (D+  ♦  D_)E  f  and  FQ  . 
t  e  [1,«] 

The  1st  statement  follows  directly  from  estimates  given  in  the  papers  cited  above.  We 
see  that  (2.2)  is  an  additional  boundary  condition  at  t  ■*  *  of  the  rank  r+  +  (rQ  -  rQ). 

Now  we  investigate  the  decay  properties  of  Up  in  dependence  of  the  decay  properties 
of  f : 

Theroem  2.2.  If  f  fulfills  (2.12)  then  the  following  estimates  hold  for  t  >  6: 

(2.15)  I  (H  f )  ( t )  I  <  const.  ID  E_1fll4.  . 

♦  +  lt,»J 

(2.16)  I (Hgf ) ( t) I  <  const,  t  E  max  Is *a+1 ^I+eD  E  'f(s)l 

s>t 

with  e  >  0  . 

If  1  >  0  then 

(2.17)  I (H  f ) ( t ) I  <  const,  t  ^  max  Is^D  E  'f(s)l 

6  <  a  <  t 

all  constants  are  independent  of  f  and  5  . 

The  first  two  estimates  have  been  proved  by  M.  Lentini  and  H.  B.  Keller  (1980). 
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The  third  estimate  follows  from 


X  a+ 1  X  a+ 1 


■  /..>!  (  0+1  a+1  (la+1  o+1.  i-1  a-y  .  i  ,  y_  -1, 

■u  (t)l  <  c.max(e  J  e  ( t  -s  )  s  dsj  max  Is  D  E  ( s ) I 

“  i«1(1)k  S  «<s<t 


where  -X  is  the  largest  real  part  of  eigenvalues  of  J  with  X  >  0  and  k  is  the 
maximal  dimension  of  the  associated  Jordan-blocks. 

By  applying  Taylor's  formula  to 

X  a+1 

/tea+1s  (to+1  .  8“+V-V-TTds 

gi(t)  “  — - _X_  o+l - 

ea+’  t-Y 

it  is  easy  to  conclude  that  lim  g.(t)  =  const.  In  particular  Theorem  2.2  tells  us  that 

t-*»  1 

inhomogenities  which  converge  to  zero  as  a  power  function  produce  a  particular  solution 
which  converges  as  a  power  function  whose  exponent  is  increased  by  (a  +  1)r  . 

Now  we  want  to  investigate  exponentially  decreasing  inhomogenities. 

.03  -  1  -  1  (J  ~ 

If  J  D  $  ( s)E  f(s)s  ds  exists  we  can  substitute  H_  by  H  which  is  defined  by 

6  “ 

(2.18)  ( H _ f )  ( t)  -  <>(  t)  /fc  D_*_1(s)E_1f  (s)sads  . 

CO 

Now  we  prove: 

Theorem  2.3*  Let  J-  consist  of  Jordan  blocks  belonging  to  the  same  eigenvalue  -X  and 

let  k  be  the  dimension  of  the  largest  of  these  blocks* 

- 1  8/0)  a+ 1  -» 

Furthermore  let  D  E  f(t)  “  t  exp(-  — -  t  JF  (t)  with  F  e  C.  ((!,"))»  &  c  R  and 

-  a+ 1  —  -  b 

u  >  0  .  Then  for  t  >  6  : 


1 


(2.19)  I (H_f ) ( t ) I  <  const,  t  exp(-  — ~  ta+1)lF_l 


if  ReX  -  oj  >  0  . 

(2.20)  l(H_f)(t)l  <  const.  exp(-  ta+1  )t^a+1  ^+®£nt*  IF_I  ^  ^ 

if  ReX  -  o)  *  0  and  0  >  -k(a+1)  The  factor  £nt  only  appears 
if  0  =  -k(a+1 )  . 

.  ~  /  to  a+1-»  k(a+1)+0_ 

(2.21)  I ( H_f ) ( t ) I  <  const.  exp( -  t  J  •  t 

if  ReX  -  u>  =  0  and  0  <  -k(a  +  1). 


All  constants  are  independent  of  f  and  6  . 

The  proofs  work  analogously  to  Theorem  2.2.  Theorem  2.3  implys  that  exponentially 
decaying  inhomogenities  produce  inhomogenities  which  converge  with  the  same  exponential 
factor,  only  the  algebraic  factor  may  change  and  a  logarithmic  factor  may  appear.  If  H 
exists  then  it  cuts  off  the  terms  of  the  particular  solution  which  are  already  included  in 
♦(t)D_£. 

Assume  now  that  J"  consists  of  more  than  one  Jordan  block  with  different  eigenvalues 
and  that  D_E_1f(t)  has  the  form  as  in  Theorem  3.  Then  H_  and  H_  may  be  used  in  order 
to  gain  a  particular  solution  which  decays  as  fast  as  possible  according  to  the  different 
cases  of  Theorem  2.3.  Doing  this  D_  has  to  be  split  up  into  the  projections  onto  the 
direct  sums  of  the  invariant  subspaces  associated  with  different  eigenvalues  with  negative 
real  part  and  H_  resp  H  have  to  be  applied  to  the  resulting  subsystems.  We  will  call 
the  resulting  operator  H  .  Its  composition  depends  on  the  decay  properties  of  f  and 


(3.13)  u'  -  t“  J(t)u  -  t'V’ftt),  1  <  t  <  « 

(3.14)  u  €  C( [1,“] ). 

For  the  following  we  need  the  definition  of  an  asymptotic  series.  A  function  P(t)  is 


said  to  be  represented  asymptotically  by  a  formal  serieB: 

oo 

-i 


(3.15) 
if 

(3.16) 
Therefore 


P(t)  ~  l  Pt  t'1  , 


t  ♦ 


i  -  0 
.-il 


tm[P(t)  -  J  P  t  ]  ♦  0  for  t  ♦  •  and  m  >  0  . 
i  -  0 


n  .  . 

(3.17)  P(t)  «  l  P.t'%  0<t'  )  for  t  *  -  and  m  >  0 

i  -  0 

holds.  To  get  more  information  on  asymptotic  series  see  Wasow  (1965). 

Coddington  and  Levinson  (1955)  and  Wasow  (1965)  prove  an  asymptotic  representation  of 
the  fundamental  matrix  of  the  homogenous  system  (3.13)  which  we  state  in 

Theorem  3.1.  Under  the  given  assumptions  on  J(t)  the  homogenous  system  (3.13)  possesses 
a  fundamental  solution  of  the  forms 


Here  D  -  diag(d„ , • • • ,d  ),  Q(t) 
i  n 


4>< t)  -  P(t)  tDeQ(t) 


diag(g.(  t) ,  *»'.g  ( t ))  where  the  cr,  / 1) 
in  1 


are  polynomials 


in  t  of  degree  a  +  1 ,  so  that 

Q(t) 

with  diagonal  matrices  Q^.  P(t) 

P(t) 


a+’  a 

J  -  +  Q.  —  +  ••.  +  Q  t 
0  a+1  la  a 


admits  the  asymptotic  expansion: 


~  I  P.t  1  for  t  ♦  » 
i  =  0 


and  PQ  =  P(~)  =  I  . 

The  unknown  coefficients  P^,  and  D  can  be  calculated  by  algebraic  operations 
from  the  Jj's.  We  get  an  algorithm  for  the  calculation  of  these  coefficients  by  inserting 
the  asymptotic  series  for  $  and  J  into  the  differential  equation,  by  formal 
differentiation  of  $  and  by  comparing  the  coefficients  of  the  corresponding  powers  of 
t  ( see  Wasow  ( 1965 ) )  . 

For  an  arbitrary  matrix  B  let  diag(B)  denote  the  matrix  which  has  only  rero  off- 
diagonal  entries  and  the  same  main  diagonal  as  B  . 
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In  the  case  a  -  0  we  find  D  -  diagtJ^).  For  a  >  0  we  pursue  in  the  following 


1)  Set  Q.|  diag(J^)  and  let  P^  be  the  solution  matrix  of  the  equation: 

pi  Jo  -  Jo  p“i  -  Ji  -  C1 

with  diagCP^ )  -  0  .  This  equation  is  uniquely  soluble  because  the  diagonal  elements  of 
Jq  are  distinct. 

2)  Pfc  and  are  determined  recursively.  Set 


0k  -  diagt^  I  i  *k_t  -  Pk.4  Qt)  ♦  Jk)  for  2  <  k  < 


and  let  Pk  with  diag(Pk>  ”0  be  the  solution  matrix  of: 


Pk  J0  "  J0  \  ^  <Ji  \ -t  "  Pk-4Ci)  +  (Jk  *  °k) 


3)  Set 


The  matrices  JQ,  •••,  Ja+1  determine  Q0,*,,,Oa  and  D  .  In  the  presence  of  eigenvalues 
of  At®*)  with  a  real  part  zero  it  is  therefore  not  sufficient  to  know  At"*)  -  AQ  in  order 
to  determine  whether  the  fundamental  solutions  which  correspond  to  these  eigenvalues  are  in 
C( [1 ,“] )  or  not.  More  (at  most  a+2 )  coefficients  of  the  series  expansion  of  Att) 
have  to  be  known.  Moreover  an  eigenvalue  with  real  part  zero  can  produce  a  basic  solution 
which  is  exponentially  increasing  or  decreasing,  algebraically  increasing  or  decreasing, 
constant  or  undampenedly  oscillating  as  t  approaches  infinity.  So  the  i-th  basic 
solution  is  in  C([1,“*J)  if  and  only  if: 


q  (t)  d 
lim  e  t 
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Therefore  we  require  at  most  that 


(3.22) 

lf(t)l  -  0(t“a_1"e)  e  >  0  . 

Now  we  de  f ine : 

r0 

(3.23) 

(a) 

u  »  u  +  u  +  u  ,  (b)  u  -  > 

p  p+  p0  p-  po  j :  i 

u 

P0i 

(3.24) 

(a) 

“p 

”  Hf,  up  “  H+f,  Up  «  H0f,  Up  -  H_f, 

(b) 

u  * 

P0i  01 

and 

(3.25) 

u 

I 

(t)  -  4>( t >  /fc  D+$~1(s)s“E"1f(s)ds 
+  08 

(3.26) 

u  (t)  ”  $(t)  D  ^  1(s)8°E  1f(8)ds 

P-  S  ■ 

(a) 

r  ♦( t)  /fc  Doil(i_1(s)saE-1f(s)ds  ,  if 

1 

(A)  holds 

(3.27) 

u  (t) 
p0i 
(b) 

-  < 

4>(t)  /*  D  ♦“1(«)aV1f<s)as  ,  if 

6 

(B)  holds 

where : 

(A) 

Re  qVi 

(t) 

♦  “  or  (Re  q^  +i  i  0  and  >  -Re 

V* 

(B) 

Ns  qr++1(t)  4 

•  -*•  or  (Re  q  ,.50  and  e.  <  -Re 
r++i  i 

V1’  • 

Then  u^”)  exists  and  fulfills  the  linear  equation: 


-1 , 


Jn (D  +  D  )u  (»)  -  -(D  +  D  )E  f(»> 

u  +■  —  p  +  — 

We  prove  that  u  is  continuous  in  t  *  "  ,  the  proofs  for  u  and  u  are  similar. 


Using  Theorem  3.1  and  (3.28)  we  get 


u  (t)  «=  P(t)tDeQ(t)  /*  D  e”Q<8)s'D+aIp"1  (slE'Vfslds  . 
P+  -1  * 

Looking  at  the  i-th  component  of  P  (t)u  ( t)  we  find, 


(P-1(t)up  (t))1  «  /te  1  1  (^)  1  s°‘(P-1(s)E~1f(s))ids 


ft  qi(t)  ’  qi(S,rtA  a  -1.  -i 


Applying  del’Hospital' s  limit  theorem  to 
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Jo  S'-  — 


(P  (o)E  f(8))1ds 


-q^t)  -di 

B  t 


we  get  because  of  the  convergence  of  the  integral  in  the  numerator 


(p"1(t)E_1£(t)) , 


(u  (“)).  -  lim  g At)  »  lim  - - - — T 

p+  1  t~  1  -q^(t)t_“  -  dit  °  1 


(E_1f  (•»)  )  j 


We  hove  used  that  P(t)  «  I  +  0 ( t“ 1 )  and  P-1(t)  »  I  +  0(t_1)  are  continuous  in  t  -  “  . 
As  in  paragraph  2  we  substitute  uh  +  up  into  the  boundary  condition: 

(8^(1)  <D0  +  D_)  +  Bo>ED0Q)5  -  6  -  B1EU  ( 1 )  -  BJfti  («)  . 

There  is  the  projection  onto  the  direct  sum  of  invariant  subspaces  associated  with 

those  imaginary  eigenvalues  of  Jq  ,  which  produce  basic  solutions  which  do  not  vanish  at 
infinity. 

Theorem  3.2.  The  problem  (3.1),  (3.2),  (3.3)  under  the  assumptions  (3.4),  (3.5),  (3.10) 


Vr- 

has  a  unique  solution  y  for  all  f  fulfilling  (3.21)  and  S  e  R  if  and  only  if 

rank[B1E*(1)(D()  +  D_)  +  B^ED^]  -  rQ  +  r_ 

where  B.  and  B  are  (  r  +  r  )  x  n  matrices.  This  solution  depends  continuously  on 
1  00  0 

*  -1  -1 

B.  (D+  +  D_)E  P  f  and  FQ  where  Fq  fulfills  DoiFg  =  FQi. 

For  the  proof  of  the  continuity  statement  we  need 
Theorem  3.3.  If  f  fulfills  (3.21)  then  the  following  estimates  hold  for  t  >  6: 

(3.28)  l(H+f)(t)l  <  const.  ID^'1  E_1  f I  ( ^  . 

If  y  >  0  then: 

(3.29)  l(H_f)(t)l  <  const.  t~T  max  ls1fD_P-1E_1  f ( s) I  . 

5<s<t 

For  i  -  1(1 )rQ: 

If  Re  q  .  .  ♦  then 

r+  +  1 

(3.30)  I(H0if)(t)l  <  const.  tt  • 
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marntMilmMk 


If  Re  qf  +  ^  ■  0  and  >  -Ra  df  +  than 


(3.31) 

l(H0if)(t)l  <  const,  t  ,rll(t(t]  . 

If 

Re 

♦  -*•  and 

Y  >  0  then 

(3.32) 

l(Hfl.f)(t)l  <  const,  t  Y  max  lsYF.(s)l  . 

«<s<t 

If 

Re 

V  +  i 

3  0  and 

e,  <  -Re  d  .  then 

i  r  +  i 

(3-33)  l(H01f)(t)l  <  conat.  t  IF^^  fcJ 

If  Ra  g  .  .  =  0  and  e,  m  -Re  d  ,  then 
r+  +  i  i  r+  +  i 

-ei 

(3.34)  l(HQif)(t)l  <  const,  t  fntlFj^ljj  . 

All  Constanta  are  Independent  of  6  and  f  . 

He  prove  estimate  (3.30)  and  the  others  are  proven  quite  analogously. 
From  (3.12)  and  (3.27)  we  derive 


f°  exp(-Re  ar  +  ^s))* 


-Re  d  ^  a  -  k. 
r  +  i  1 


ds 


IH0if,(t).  <  const. - 


-Re  d 


exp ( -Re  qr  +  1<t))t 


r+  +  1 


'Fi'[t,-) 


0+1  -k^ 

Applying  del1 Hospital' s  rule  and  using  |Re  qr  +  A ( t > |  “  o(t  )  the  estimate  (3.30) 

follows.  Herewith  Theorem  (3.2)  is  proven  completely.  An  inhomogenity  f  fulfilling 
(3.22)  produces  a  particular  solution  which  decays  at  least  as  fast  as  t_E  int. 

Now  we  again  investigate  exponentially  decaying  inhomogenities . 

We  define  H  as  in  Chapter  2  ,  H  as  (3.27a),  H  as  (3.27b)  no  matter  whether 
(A)  or  (B)  holds.  Of  course  the  definitions  make  sense  only  *or  inhomogenities  which  make 
the  appearing  integrals  exist. 

We  assime  that 

(3.3?)  f(t)  -  ep(t)tSF(t),  F  e  cytl,-)),  6  e  R 
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where  p(t)  is  a  real  polynomial  of  maximal  degree  a  +  1  and  p(t)  +  -*  for  t  ♦  “. 
We  split  up  the  operators  H_,  H_  into  H_j  ,  H  by  substituting  the  projection  D_ 

by  D_j  (for  j  -  1(1)r_)  which  is  the  projection  onto  the  invariant  subspace  generated 
by  the  j-th  eigenvalue  of  Jg  with  real  part  smaller  than  zero. 


Theorem  3.4.  Let  (3.35)  hold.  Then  the  following  estimates  can  be  derived  for  t  >  6: 


8+a+l-t. 

I(U  .  f)(t)l  <  const.  epl  't  IFI., 

-3  [o/t] 


if  p(t)  -  Re(  c 


r+  +  r0  +  j 
the  polynomial  p(t)  -  Ret 


.  <  t) ) 


and  1 ^  is  the  degree  of 


\  +  r0  +  j 


(3.37)  l(H  f)(t)l  5  const.  epl t! ta+8+1 »F» . . 

-3  [6,t] 

if  p(t)  -Ref  +  ^  +  .  ( t ) ) 


and  8  >  Re( d  )-a-1 . 

r  +  r„  +  i 
+  0  J 


(3.38) 


l(H  ,f)(t)l  5  const.  eP^t^ta+®+1  £nt  IFI 


[6.t] 


if  p(t)  -  Re(  q^  ___  ^  ^  ( t ) )  and  6  =  Re(d_  ,  _  _  ^)-a-1 


r>  +  ro  +  j 


r+  +  r0  +  j 


(3.39)  l(H  f)(t)<  <  const.  eP( t>ta+8+1 ,Fl 

-3  [t,'»] 


if  p(t) 


Re(c 


( t ) )  and  6  <  Retd  .  >-a-1. 

ro  +  j  rt  +  r0  +  ^ 


_  B+a+l-f, 

(3.40)  I  (H  f )  ( t )  I  <  const.  eP'  't  L|FI, 

-j  It,“J 

if  p(t)  -  Re(q  .  (t))  ♦  . 

+  0  J  ^ 

The  estimates  (3.37),  (3.38),  (3.39)  are  also  valid  if  resp.  H  and  the  index 

r+  +  rQ  ♦  j  are  substituted  by  H__.  resp  and  the  index  r+  +  i  for  i  =  3(1)  Tq  . 

The  proof  of  Theorem  3.4  is  quite  similar  to  that  of  Theorem  3.3. 

As  in  the  constant  coefficient  case  exponentially  decaying  inhomogenities  produce 
particular  solutions  which  convorge  with  the  same  exponential  factor,  only  the  algebraic 
factor  may  change  and  Int  as  a  factor  may  appear.  If  the  inhomogenity  f  contains 
(int)s  with  s  c  N  as  a  factor  then  all  derived  estimates  are  valid  if  their  right  hand 
sides  are  multiplied  by  (fnt)3.  This  follows  by  direct  estimation  of  the  appearing 


integrals  and  by  using  the  Theorems  3.3  and  3.4.  The  operator  which  produces  the  best 


order  of  convergence  according  to  Theorems  3.3  and  3.4  and  which  is  composed  of  the 

operators  H  ,  H  H  ,  is  again  denoted  by  H.  The  composition  of  H  depends 

+  01  0i  -i  -i 

on  the  decay  properties  of  f  and  on  JQ,  J^ ,  •*•  ,  • 

We  want  to  weaken  the  condition  (3.6)  which  is  a  very  strong  assumption  on  Aft).  The 

basic  idea  for  this  is  that  the  determination  of  the  matrices  ,  D  only 

requires  the  knowledge  of  A  .  •••  ,  A  .  while  the  matrices  A  ,  •••  do  not  influence 

0  <X+ 1  Ct+2 

the  exponential  and  algebraic  factors  of  the  fundamental  matrix,  they  only  influence  the 
power  series  factor  P(t). 

We  assume  that 


(3.41) 

and  therefore  we  can  write 


A(  i  )  «  C“+2(I0,  ’))  , 

A  e  C(  [1,®]  ) 


6  >  1 


(3.42)  Aft)  »  ft0+t_1A  +  t”®"1*  +1  +  A(t),A(t)  =  a(t)t 

where  a  e  C.  ((1,“))  and  8  >  0  . 

D 


The  equation 


y'  «  taA(t)y  +  f(t)  ,  t  >  1  ,  ye  c([l,»)) 

where  f  fulfills  (3.22)  can  be  rewritten  as 

01+1  i 

(3.43)  v*  «  t“(  l  A  .t  )y  +  (A(t)y  +  f  ( t) )  ,  y  e  C({1,»J). 

1  »  0 

The  homogenous  problem  (3.43)  with  f  i  0  has  the  general  solution 

(3.44)  yh  “  E  ♦<®o  +  D.,?  +  E  ”l*yh'  5  e  C" 

where  E  transforms  Ag  into  its  Jordan  form  and  E$  is  the  fundamental  matrix  of  the 
unperturbed  problem 


y 


h 


a  +1 

t°(  [ 


Now  let 

(3.45)  l$(  t)  (D  ♦  D  )l  -  p(t)tdea(t)  -  p(  t)a  ( t) ,  p  e  C([1,«]) 

O  “  n 

H ^  Is  composed  so  that  inhomogeni ties  which  decay  as  t  a  2  ®o^(t)  produce  fast  decaying 
particular  solutions  regarding  Theorems  3.3  and  3.4. 
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From  (3.44)  we  get  the  equation: 


(3.46) 

for  which  we  take  as  basic  Banach-space: 

(3.47) 


(I  -  EH  A ) y  »  E$(D  +  D  >? 
l  h  u  - 


(A„  ,  -  { u  t  u  ( t )  -  U(t)o.  tt),  U  E  C.  <[«,»))},  lul  ,  -  lul „,). 

0^,6  h  b  0^,6  [o,*] 


If  o  (t)  S  1  then  we  set  A  *  »C([fi,0#)). 

h  ch.5 

we  get  the  estimate 

(3.48)  IH.AI  ,  •  max  *H.Ay.  I  ,  <  const.  6  1  ®Jtn6  <  - ;  ~ 

1  °h'5  ly  I  ,<1  1  h  V*  2,E‘ 

n  o  ,o 
h 

if  i  is  sufficiently  large.  Therefore  (I  -  EH ,A )  1  exists  as  an  operator  on  A  . 

1  V5 

and  for  £  e  C  : 


(3.49) 


yh  -  (I  -  EH1A)'1E$(D0  +  Djc  -  ^(Dfl  +  DJ5  £  Afl  &  . 

h' 


As  particular  solution  yp  of  (3.43)  we  set. 


(3.50) 

The  inhomogenity  f  fulfills 

(3.51) 


yp  -  EH3Ayp  *  EH^ 


lf(t)l  ■  0(^0^^*)  and  op(t)  -  ta+1+dint  eq*t^  ♦  n 


and  H2  is  composed  ao  that  (HjfHt)  decays  as  fast  as  possible  with  regard  to 


Theorems  3.3  and  3.4.  Then 
(3.52) 


I (H  f ) (t) ■  -  0 <0  <t)) 
2  P 


is  composed  to  make  particular  solutions  belonging  to  inhomogenities  which  decay  as 
t  a  2  (t)  decrease  as  fast  as  possible.  As  basic  space  we  now  take  A  .  and 

p  _  __  Op-* 

conclude  the  invertibility  of  (I  -  EH  A)  on  A 


with  6  sufficiently  large  and  get 
(3.53) 

Obviously 


0  ,6 
P 


y  -  <ji(f)  -  (I  "  EH,A)'1EH,f  E  A„  . 
p  3  2  0  , 0 

P 


(3.54) 


y(t)  »  y.  f t: )  +  y  (t)  c 
h  p 
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holds.  By  substituting  H  ,  which  is  defined  by  (3.24),  for  in  (3.44)  it  is  easily 

shown  that  the  solution  manifold  yh  (with  the  parameters  (Dq  +  D  )5 )  is  unique  in 

C< [4,«] )  ,  because  A  .  t  C([6,»])  and  because  the  solution  space  is 
0.  ,0 
n 

rQ  +  r  dimensional.  Therefore  y  defined  in  (3.54)  is  unique  in  C(  [  6  ,«•)  )  (as 
manifold) . 

In  order  to  get  the  solution  in  C([1,*l)  we  solve  the  'regular'  problemi 

(3.55)  y'  -  taA(t)y  +  t®f(t)  1  <  t  <  4 

(3.56)  y<4)  -  y.  (4)  +  y  (4)  . 

h  p 

From (3. 44)  and  (3.49)  we  get,  using  the  expansion 

00 

(I  -  G)"1  -  l  GL  for  IGI  <  1 
i  -  0 

the  following  estimates  which  hold  for  t  >  4  s 

(3.57)  li|i°(t)  -  E$(t)  (D.  +  D  )  I  <  const.  t“1-B  tnt  •  a.(t) 

-  0  -  n 

and 

(3.58)  I  (*(  f ) )  ( t)  -  E(H  f)(t)l  <  const.  t“1_B  tnt  •  o  (t>  . 

2  p 

Theorem  3.2  remains  valid  if  E$(1)(D0  +  D_ )  is  substituted  by  i)i°(1)  where  i)i®(t)  has 
bean  continued  to  [1,**]  . 
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4.  Linear  Variable  Coefficient  Problems  -  General  case. 

In  this  paragraph  problems  of  the  form  (3.1),  (3.2),  (3.3)  fulfilling  (3.4),  (3.5), 
(3.6)  but  not  (3.10)  are  investigated.  So  we  assume  a  general  Jordan- structure  of  Jg  in 
the  block  diagonal  form 

(4.1)  JQ  -  diag(jJ,J°,j"),  dim(jJ)  -  r+,  dim(J°)  =  rQ,  dim(j")  =  r 

where  J*,Jg  resP-  Jg  contain  the  Jordan  blocks  which  are  associated  with  eigenvalues  with 
real  part  greater,  equal  reap,  smaller  than  zero.  As  in  the  case  of  distinct  eigenvalues 
of  Jq  we  determine  an  asymptotic  expression  for  the  fundamental  matrix.  This  expression 
is  given  by  Wasow  (1965): 

Theorem  4.1.  Under  the  given  assumption  on  J(t)  there  is  a  fundamental  matrix  which  has 
the  form 


4><t) 


P(t)t 


D  Q(t) 
e 


where  Q(t)  is  a  diagonal  matrix: 


Q(t)  *  diag(J  ) 


a+1 


a+1- 


a+1  -  — 

-  P 


0'  a+1  ’  V1  p( a+1 )-1  +  ®2  p(a+1 )-2 

2 

P  — 

•  •  •  +  o  — — —  +  o 

vp(a+1)-2  2  p(a+1 )-1 


with  some  p  c  N  ,  ,  D  is  a  constant  matrix  in  Jordan  canonical  form  and 

P(t)  =  P1 (t)P2(t) 

where 

» 

P^t)  -  I  ♦  ^  l  i  Put_i  ,  T  +  - 

and 

-  _  1 

P2(t)  ~  l  P  t  P,  t  +  - 
i  =  0 

The  (diagonal)  elements  of  Q(t)  which  correspond  to  a  particular  Jordan  block  of  D  are 
equal.  Therefore  tD  and  e®*fc'  commute.  Moreover  the  block  structure  of  D  is  a 
subdivision  of  that  blocking  of  Jg  which  is  obtained  by  aathering  all  Jordan  blocks  of 
J0  belonging  to  the  same  eigenvalue.  Also  Pjft)  has  a  block  structure  which  is 
identical  to  the  above  mentioned  blocking  of  Jg  . 


-20- 


has  the  following  form 


k-1 

(4.5)  p  (t)  =  t  (t)  n 

i*1 


ri  +  +  ri 


Ti+i(t) 


,  p,(->  -  i 


Because  of  (4.3)  k  systems,  where  the  leading  matrix  of  each  of  them  has  only  one 
eigenvalue,  have  to  be  treated  now. 

2)  We  substitute 

a+1 

u(i)  -  v(i)  expK  ) 

and  get 


(4.7) 


V  * 

(i) 


ta(J(i)(t) 


Vr^i) 


for  i 


The  leading  matrices  of  the  systems  (4.7)  are  now  -  Hi 


3)  We  apply  so  called  shearing  transformations 
(4.8)  v 


'(i) 


S(i)(t)w(i) 


=  1(1)k  . 

having  the  only  eigenvalue 


where 

-g  ,  -2g  -  ( r  - 1  )  g 

(4.9)  S(i)(t)  *  dla9(l.t  t  , 


with  g^  rational  and  g^  >  0  to  the  systems  (4.7).  The  g^  should  be  chosen  so  that 
the  leading  matrices  of  the  resulting  systems,  which  have  w(ij  as  dependent  variables, 
have  more  than  one  different  eigenvalue  or  if  not  possible  for  a  certain  i  ,  so  that  the 
rank  of  this  new  i-th  system  is  smaller  than  a+1  or  that  this  systems  splits  up  into 
separate  subsystems. 

Wasow  (1965)  showed  that  it  is  always  possible  to  achieve  one  of  these 
simplifications.  In  order  to  get  systems  where  only  integral  powers  of  the  independent 
variable  occur  we  substitute 

1/(g-a-1)  1/p 

(4.10)  Xi  =  pi  '  t  for  i  *  1 ( 1 )k 

where  p^  is  the  smallest  Integer  so  that  g^p^  Is  a  natural  number.  Then  we  get  systems 
of  the  form 


(4.11) 


'<W 


xi  C(i!,xi,w(i),Xi) 


..  ■  ..  •  I'xiiUBiinA’ii  i 


where 


(4.12) 


(a)  hL  -  (a  +  1  -  g.Jp^l,  (b)  CU)‘V~  ^  C(i)Xi 


-3 


4a)  If  has  only  0  as  eigenvalue  then  we  have  reduced  the  rank  of  the  system  or  it 

splits  up  into  separate  subsystems  of  lower  order.  Applying  more  shearing  transformations 
to  (4.11)  we  end  up  with  a  system  whose  leading  matrix  has  either  more  than  one  different 
eigenvalue  or  has  rank  equal  to  0 .  In  the  second  case  we  have  a  system  with  a  singularity 
of  the  first  kind  for  which  the  fundamental  matrix  $  (t)  has  the  form 

U.13)  ■  p(i,<t>tD<i>'  p(i)(t)  '  J  /a)*'* 

and  is  a  constant  matrix  in  Jordan  canonical  form  (see  WaBow  (1965)). 

0 

4b)  Now  we  assume  that  C  ^  has  at  least  two  different  eigenvalues.  Then  we  transform 

0  ~  0 
to  its  Jordan  canonical  form 


(4.14) 

and  substitute 

(4.15) 


0  ~0  — 1 
C(i)  "  E(i)C(i)  E(i> 


w(1)(Xi)  -  E(i) 


~0 


getting  a  system  whose  leading  matrix  is  C, ,  that  means  in  Jordan  form 

hi 

<4.16)  2;i,(v  -  xi  3(i)<xi)2(i,(xi)  • 

5)  We  apply  the  transformation  explained  in  1)  to  the  system  (4.16)  in  order  to  get 

separate  subsystems  of  lower  order  whose  leading  matrices  have  the  only 

eigenvalue  •  By  bhe  means  of  2)  we  normalize  these  systems  so  that  their  leading 

3 

matrices  have  only  the  eigenvalue  0  .  These  transformations  split  off  the  factors 


V1 


(4.17) 


exp^(i)j-hpr  )  for  i’1(m 


Resubstituting  Ln  (4.10)  and  using  (4.12}(a)  we  notice  that  the  argument  in  (4.17)  is  of 
a+1 -g^  a+1 

the  order  t  that  means  of  order  lower  than  t  which  is  the  order  of  the 

argument  of  the  first  exponential  factor  because  if  g^  »  0  the  system  remained  unchanged. 


Applying  another  set  of  shearing  transformations  3)  w.  arrive  at  4a)  or  4b). 
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6)  So  a  finite  chain  of  all  the  in  1 )  -5 )  described  transformations  result  in  a  set  of 


one  dimensional  systems  and  systems  with  a  singularity  of  the  first  kind.  Setting 
-m  - 

p  -  p  where  p  in  the  smallest  common  multiple  of  all  the  p's  used  in  the  sequence  of 
shearing  transformations,  and  m  is  the  number  of  these  transformations  which  split  the 
system  into  a  set  of  systems  described  above,  we  get  the  formula  for  the  fundamental  matrix 
$(t)  given  in  Theorem  4.1  by  taking  into  account  (4.13).  Moreover  we  get: 


(4. 18) 


P2(t) 


(.n.Si(t,EiP2i(t))P4(t)'  P2i(t) 


I  + 


i=1 


l  « 

j  =  1 


.  1 
-  P 


2i 


and  the  S^tt)  are  composed  of  submatrices  S^^ft)  defined  in  (4.9).  The  E^’s  are 
regular  and  P4(t)  is  derived  by  solving  the  systems  with  singularities  of  the  first  kind 
using  (4.13). 


We  define: 

m 

(4.19)  P3(t)  =  n  Si(t)EiP2i(t) . 

-1  1  =  1 

An  estimate  of  p^  { t )  can  be  obtained  in  the  following  way. 

Let  be  the  projection  onto  the  direct  sum  of  eigenspaces  associated  with  the 

eigenvalue  X^  of  JQ  .  Then 


(4.20) 


IP2^(t)D^I  <  const,  t 


[(ri'1)(gii+  77  * 


9i2  .  gi3 


im 


pll”'pim-1 


)] 


holds.  The  sum  in  the  exponent  of  t  is  derived  by  estimating  s”1 (t)  and 

by  taking  into  account  the  block  structures  of  the  E^1  and  P2^(t)  .  p^ 

and  g^j  are  as  in  (4.10)  and  represent  that  sequence  of  shearings  starting  off  from  the 

i-th  r .-dimensional  subsystem  and  givina  the  largest  exponent  in  (4.20). 


For  this  sequence  we  can  calculate  the  ranks  of  the  corresponding  seouence  of 


subsystems  as  in  (4.12) 


(4.21) 


hi0 

+ 

1 

«  a  +  1 

hi1 

+ 

1 

-  P^Chio  *  1 

-  gt1) 

hi2 

+ 

1 

’  Pi2(hi1  +  1 

’  gi2’ 

h  +  1 
i  ,m 


Pin/hi,m-1  + 


gim> 


By  assumption  h,  +  1  >  0  holds  and  we  get  recursively: 
1 ,  m 


(4.22) 


yi2 

a  +  1  >  g,  +  -  +  •  •  •  +  — 

11  Pi1  P 


lm 


il'  pi,m-1 


and  therefore  the  estimate 

( r  “ 1 ) ( o+ 1  ) 

IP3  ( t)D^ I  <  const,  t 

holds. 

The  basic  solutions  with  4>{  t)  -  ( t) ,  •  •  •  ,V>1  ( t) )  fulfill 

Vt)  d  j 

(4.24)  l^i(t)«  <  Pl(t)e  1  t  Unt)  ,  pt  e  C(  [1,-]  ).  • 

Eigenvalues  of  Jg  with  positive  real  part  produce  exponentially  increasing,  eigenvalues 
with  negative  real  part  produce  exponentially  decreasing  basic  solutions.  Imaginary 
eigenvalues  of  Jg  can  produce  exponentially  and  algebraically  increasing  and  decreasing 
constant  and  oscillating  and  logarithmically  increasing  solutions.  The  asymptotic 
behaviour  of  a  particular  basic  solution  can  be  determined  by  the  knowledge  of  Q(t), 

D  and  P  ,  •••,?  where  m,  is  sufficiently  large.  Therefore  the  solution  of  the 
homogenous  problem  (3.16),  (3.17)  (under  the  assumptions  of  Chapter  4)  is: 

(4.25)  u^t)  -  $(t)(Dg  +  D_K  ,  C  E  c" 
where  DQ  and  D_  are  defined  as  in  Chapter  3. 


Nov  we  construct  a  particular  solution  of  the  problem 


(4.26) 

(4.27) 

We  substitute 


u'  «  t°J(t)u  +  t°E  1f(t)  1  <  t  <  K,  f  e  C([1,“]) 

P  P 

u  e  c( [1,«] )  . 

P 


(4.28) 


and  define 


(4.29) 


u  (t)  -  P„(t)v  <t)  and  v  (t) 

pip  p 


v  (t) 


V  (t) 
P0 

V  (t) 
p_ 


V  (t) 

p+ 

0 

u  (t)  -  P,(t) 

p+  1 

0 

,  u  (t)  «  P  (t) 

po  1 

v  (t) 

po 

0 

r- 

0 

u  (t)  «  P  (t) 
p_  ' 


Vp  (t) 


so  that  the  definitions  (3.23a)  and  (3.24a)  hold.  We  get  the  separated  problems 


— 

• 

r  v  i 

p> 

V 

po 

-  t“(diag(J*  +  J+(t),J°  +  J°(t),  j‘  +  J_(t)) 

< 

►0 

o 

V 

V 

P_ 

L  P-J 

+  t°P^1(t)E*1f(t) 


where  J+(t),  J®(t),  J-(t)  have  an  asymptotic  power  series  expansion  in 

t-1  without  a  constant  term.  We  define  the  block  components  v  ,v  as  in  de  Hoog  and 

p+  p- 

Weiss  ( 1980a, b)  as  solutions  of  the  operator  equations: 

(4.30 )( a)  v  =  H  J+v  +  H  (P? 1 f ) 

P.  ♦  P*  +  1  + 


(4. 30 ) ( b) 


*  -1  -1 

V  =  H  J  v  +  H  (P,  E  f) 

p  -  p  -  1 
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where  (P^E  1f)+  resp.  (P^E  1f)_  are  the  first  r+  reap,  the  last  r_  components 


of  P^E”1f  and  H+,  H  are  the  operators  defined  similar  to  (2.11). 


(4  •  30  )  ( c) 


(H+g)  ( t)  -  /texp(^1  (t°+1  -  sa+1))sag(s)ds 


(4.30 )( d) 

for  €  C(  [1,*>]  )■ 


j 

*  ....  rt  r  0  ,  a+1  a+t.'v  a  .  .  , 

(H_9)(t)  »  J  exp[  —  (t  -  s  )Js  g(s)ds 

6 


From  (4.30)  we  derive 

(4 .31 ) (a)  v  -  (I  -  H+J+)“1H+(p”1E-1f)+  E  C([5,®]) 


(4.31 )  (b) 


-,-1„  ,.-1-1. 


v  -  (I  -  H  J  )  H  (P  E  f )  E  C([6,»]) 

P_  -  -  1 

with  S  sufficiently  large.  The  proof  of  the  invertibility  of  (I-Hf,-J+,“)  is  given  in 
de  Hoog  and  Weiss  (1980a, b).  Again  we  get: 


(4.32)  (a)  v  (-)  -  (J*)"1(E-1f(-))^,  (b)  v  (»)  *  (J j'1(E_1f(»)) 

P+  0  +  p_  0 

Jbecause  J+(“)  -  0  and  J-(»)  -  0  . 


DQP1"1(t)E'1f(t)  -  F0(t:t"r(a+1)'e  ,  FQ  e  Cb([1,»)) 


Now  we  assume  that  for  some  e  >  0 : 

(4.33) 

holds,  where  r  is  the  largest  algebraic  multiplicity  of  the  eigenvalues  of  JQ  with  real 
part  zero. 

So  r  is  defined  differently  to  r  in  Chapter  2.  However,  more  sophisticated 
assumptions  on  f  could  be  used  (similar  to  (3.21))  but  we  will  use  (4.33)  for  simplicity. 
The  system 

(4.34)  v'  -  t“(J°  +  J°(t))v 

Pn  0  Pn 


is  composed  of  separate  systems,  each  of  them  associated  with  one  imaginary  eigenvalue  of 
Jq  and  r  is  the  maximal  dimension  of  these  subsystems. 

We  take  one  of  these  (inhomogeneous)  subsystems 

vp0(i)  '  taj(i)!t,vPo(i)  +  ta(p;1(t)E-1f(t))0(i) 


(4.35) 


*  -1  -1 

where  -  (P  (t)E  f  { t)  >  0  (  i )  consists  of  the  corresponding  components  of  the 

inhomogenity  P-1  (t)E-1f  (t) .  The  leading  matrix  J*  of  J^j(t)  has  the  only  eigenvalue 
X^  with  ReX^  -  0  • 

Now  we  apply  the  transformations 


(4.36) 


t0(i) 


X  t“+1  “ 

=  )  *  S(i)<t)WpQ(i)  '  Xi=  CU)tPi 


(4.37) 


w  (x  )  =  E  z  (x.  ) 

P0(i)  1  (i)  P0(i)  1 


as  defined  in  (4.6),  (4.9),  (4.10)  and  (4.15)  and  get 

h 

(4.36)  z 


i~  *V 

(x  )  »  x.V  ..  (x,  )z  (x  )  +  x  g,..(x.) 

p0(i)  1  1  (i>  1  P0(i)  1  1(1,1 


where 


-  *>  Pi  9.  .  .  x  p  X  x  p.  (o+l)  .  x.  p 

(4.39)g(1)(xi,  -  (Suj-J  ■u)<i)((^ir)  )‘W^(TT)  ) 


,g(i)<Xi)'  "  const  e 


-p.e-rpi(a+1-gi) 


From  (4.23),  (4.33)  we  derive: 

(4.40) 

If  the  leading  matrix  of  the  system  (4.38)  has  at  least  one  eigenvalue  different  from  zero 

then  the  separating  transformation 

z  =  P,..,(x  )z  ,  P  (“)  »  I 

P0(i)  2(i)  1  P0(i)  2(i> 

can  be  applied.  All  resulting  sybsystems,  whose  leading  matrices  have  eigenvalues  with 
real  part  different  from  zero  can  now  be  solved  by  the  means  of  (4.30),  (4.31)  because  the 
new  inhomogenity  has  the  form 

(4.41) 

For  all  other  subsystems  this  sequence  of  substitutions  is  repeated  as  long  as  we  arrive  at 
systems  whose  leading  matrices  have  eigenvalues  with  real  part  different  from  zero  or  one¬ 
dimensional  systems  or  we  arrive  at  systems  with  a  singularity  of  rank  zero.  In  the  second 
and  third  cases  only  eigenvalues  with  real  part  zero  have  been  split  off  therefore  the 
inhomogenities  do  not  contain  exponentially  increasing  or  decreasing  factors  (see  (4.39)) 


9(i)(xi)  *  P2(i)  <xi,9(i)  <Xi>  E  * 
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and  therefore  every  inhomogenity  gg  which  occurs  in  this  sequence  of  transformations 
fulfills  after  resubstitution  to  t  as  independent  variable: 


(4.42) 


•  gQ(t)l  <  const.  IPj1  (tJDjt0*1  ,Dgp11  (t)E  1  f ( t > I  <  const,  t  E 


because  (4,22)  holds. 

A  particular  solution  for  one- dimensional  systems  can  be  found  according  to  Chapter  3, 
so  we  just  have  to  treat  systems  with  a  singularity  of  rank  zero. 


We  want  to  solve 


(4.43) 


z'  «  (1  b  +  —  B( x) ) z  +  —  g  (x),  B(x)  =  B(x)  - 
p  X  x  p  x  X 


(4.44)  z  (x)  >  z  (x)x  E(inx)^,  z  e  C.  ([1,**)>,  j  e  Nn 

p  p  p  D  0 


where  B  is  a  constant  matrix  in  Jordan  form,  B  t  C([1,“])  and  gQ  fulfulls  (4.42). 


(4.45) 


B  -  diag(B  , • • • ,B  )  ,  B 

i  q  i 


b  '-I  I  ,  i  -  1(1 )q 


and  let  D1  be  the  projection  onto  the  invariant  subspace  associated  with  bi. 


Then  we  define 


(4.47) (a) 


(4.47 )( b) 


(Gig0)(x) 


z  -  GBz  +  Gg„  ,  G“G,  +  •••  ♦  G 
P  P  0  1  q 


B  fx  i  -B-X  ,  , „  ^  , 

x  J  D  s  gn(s)ds,  b  <  -e 

«  0  1 


B  rm  i  -B-I 

x  J  D  s  g^  (  s )  ds  ,  b^  >  -e 
x 


It  is  easily  checked  that  for  x  >  6 

(4.48)  KGg  )  ( x)  I  4  const.  (lnx)^max(  max  Ig  (s)l,  x  £  •  max  ls£g  (s)») 

SE[X,»]  sc [4 ,x) 

holds  where  j  -  max  (dim  (B  )).  Therefore  we  get  from  (4.46) 

1  ~ 

(4.49)  z  -  (I  -  GB)  Gg. 

P  0 

Zp  fulfills  (4.44)  because  of  a  contraction  argument  similar  to  that  applied  ir.  (3.47), 
(3.48),  (3.49).  If  90<x)  "  0(x  C(lnx)i)  then  the  right  hand  side  of  (<*.48)  has  an 
additional  factor  (inx)1. 


1 


After  having  performed  all  necessary  resubstitutions  we  get  a  solution  v  ( t) 

p0 


fulfilling  v  («)  -  0  .  Defining  r  ,  D  in  the  same  way  as  we  did  in  Chapter  3  we  get 
P0  0  00 

Theorem  4.2!  The  problem  (3.1),  (3.2),  (3.3!  under  the  assumption  (3.4),  (3.5),  (4.1)  has 

vr- 

a  unique  solution  y  for  all  f  fulfilling  (3.21) (a),  (4.33)  and  for  all  B  e  R  if 


and  only  if 


ranklB^OMDg  +  D_)  +  B^ED^l  -  rQ  +  r_ 

where  B1  and  Bw  are  ( rQ  +  r_)  *  n  matrices.  This  solution  y  depends  continuously  on 
0  ,  (D+  +  D_)E-1P1-1f  and  FQ  which  is  defined  in  (4.33). 

The  proof  is  complete  if  the  continuity  statement  is  proven. 

Theorem  4.3.  If  f  e  C([1,«J)  and  (4.33)  holds  then  the  following  estimates  hold  for 
t  >  6: 

(4.50)  I (H  f ) ( t) I  <  const.  ID  p"1E-1f I fk  , 

+  +  1  [t,»] 


(4.51) 


l(H_f)(t)l  <  const,  t  ^ • 


max  ls^D_P  ^  ( s)E  ^(s)!  for  y  >  0. 

s  e  [5,t] 


(4.52) 


KH0f)(t)l 


j0 

<  const,  (int)  •  max(t 


se[6,t) 


.  e+(u+1 )r  -1  - 

Is  D()P1  ( s)E 


f(s)l. 


max 
s  >  t 


ls<a+1)rD0P^1(s)E"1f(s)l) 


where  j  «  max  (dim(D  ))  and  D  are  the  Jordan  blocks  of  the  matrix 
°  Ak  k  1,1 

4.1  for  which  the  corresponding  polynomials  q  (t)  fulfill  Req  (t)  S  0 


constants  are  independent  of  6  and  f  . 

Proof.  From  (4.31)  we  conclude 

00  ^ 

(4. 53 ) ( a)  v  (t)  -  l  ((H+J+)iH+(P~1E_1f)+)(t) 

p+  i  «  0 

00 

(4.53Mb)  v  (t)  -  l  ((H  j')^  !  P*1  E~^  f )  )(t) 

p-  i  «  0  1 

for  6  sufficiently  large. 


D 


in  Theorem 
All 
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H  reap  H  fulfill  the  estimates  (2.15)  reap.  (2.17)  and  therefore  we  get: 

00 

(4.54)  (a)  Iv  ( t )  1  <  const.  max  ID  P  'ls)E  1  £(  s)  I  £  (■j-)* 

pt  se[t,“!  i  -  0  Z 

00 

(4.54)  (b)  Iv  (t)l  <  const.  t_1f  max  lsYD  P  -1  ( s)E_1f  ( s)  I  •  l  4)1 

P-  se[6,t]  i  -  0  Z 


if  6  is  so  large  that  < *  ( 5  “>]  *  an<J  *[6  «]  *  %  ♦  The  estimates  (4.50),  (4.51) 

follow  by  using  (4.29)  and  (4.52)  follows  from  the  derivation  of  v  and  from  (4.48), 

P0 

(4.49). 

Theorem  4.4  implies  that  an  inhomogenity  f  fulfilling 

(4.55)  f(t)  -  t'1“+1)r-£(4nt)iF(t),  Fe  C.([1,“)) 

b 

produces  a  particular  solution  Hf  for  which  the  estimate 

5  +i 

(4.56)  I ( Hf ) ( t ) I  <  const.  t-e(2nt)  °  IFI  ^ 
holds. 

Now  we  assume  inhomogen i tie a  of  the  form 

(4.57)  f(t)  -  ep<t)t6(«nt)Vt),  Fecytl,-)) 

± 

and  eP^t't®F(t)  fulfills  (4.33)  and  p(t)  is  a  polynomial  in  tP  with  leading  term 


.  It  is  possible  to  construct  a  particular  solution  Hf  “  H+f  +  HQf  +  H_f 


fulfilling 

-  j  +f 

(4.58)  (a)  ■  (HQ  f  ( t)  I  <  const,  e  p<  t)tB+rt“+1  ’  (tnt)  °  IFI[J(<>(J 


(4.58)  (b)  l(H_f)(t)l  <  const.  eP<  t!  t8+^<a+1 '  Unt)  ‘  IFI^.j 

where  k  is  the  maximal  algebraic  multiplicity  of  the  eigenvalues  of  JQ  with  negative 
real  part  and  j_  is  the  maximal  dimension  of  Jordan  blocks  in  DD_  .  To  outline  the 
proof  we  assume  the  system 
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(4.59) 


t“(J  +  J  (t)t_1)v  +  t“g(t>,  J  (»)  *  0 


where  J  has  an  asymptotic  series  expansion  and  J_  has  only  one  eigenvalue  X  with 

ReX  <  0  and  g(t)  fulfills  (4.57).  If  -ReX  +  Rep^  >  0  we  determine  a  particular 

solution  Vp  by  applying  (4.30)  (b)  .  If  -ReX  +  Rep^  <  0  we  use  the  operator  H 

instead  of  H  in  (4.30)  where  H  is  defined  by  (4.30) (a)  with  6  =  “.  In  both  cases 

P  ( t )  8 

the  estimate  Iv  (t)l  <  const,  e  'tP (int)  IF  I/,  ,  follows  by  a  contraction  argument 

P_  -  16,®] 

and  by  (2.19),  (2.22).  However,  if  ReX  =  RePg  we  perform  the  substitution  of  the  form 
(4.36),  (4.37),  (4.40)  as  long  as  we  arrive  at  systems  whose  eigenvalues  have  real  parts 
different  from  the  real  parts  of  the  corresponding  coefficients  of  p(t).  Then  we  apply 
H_  reap.  H  .  If  Rep(t)  =  Req^(t)  where  qj  is  a  polynomial  in  the  exponential  factor 
of  the  fundamental  matrix  of  the  homogenous  system  (4.59)  this  sequence  of  transformation 
either  leads  to  one  dimensional  systems  or  to  a  system  with  a  singularity  of  the  first  kind 
(rank  zero).  The  first  has  been  treated  in  Chapter  3  and  the  second  is  solved  by  (4.46), 
(4.47).  Resubstitution  and  application  to  all  separate  subsystems  corresponding  to 
eigenvalues  with  negative  real  part  yields  the  estimate  (4.58) (b).  The  estimate  (4.58) (a) 
is  gained  in  a  similar  way. 

Now  we  assume  that  the  matrix  A(t)  of  the  system  (3.1)  fulfills 
(4.60) 

A  e  C< [1,«] ) 

instead  of  (3.5),  where  l  =  max(r,k)  is  defined  for  the  Jordan  form  of  A(°°)  .  Therefore, 


(1  ’I  ( d+ 1 )  1+1  .  r .  1 .  j. 

a(;-J  c  C  (  [0,-j]  )  ,  6  >  1 


(a+1  )t 


(4.61)  A(t)  -  Aq  +  t_1A  +  •••  +  t'<a+1)SA 
where  a  e  C^( [1 ,») ) ,  B  >  0  . 

The  system  (3.1)  can  now  be  written  as : 

(4.62)  _  y*  =  t°  (A  +  •••  +  t"(a+1)^  A 

(a+1 )£ 

By  regarding  (4.62)  as  a  perturbation  of 

(a+1 )5  „ 

=  t  (  l  t  A  >y  +  t  f 
i  =  0 


+  A(t),  A(t>  -  a(t)t 


■(a+1  )f-1-B 


)y  +  ta(A(t)y(t)  +  f  ( t) ) . 


(4.63)  y 

and  proceeding  similarly  to  Chapter  3  we  find  that 


(4.64) 


liji_(t)  -  E$(t)(Dg  +  D_ )  I  <  const,  t  (Int) 


!_g  2max(j0,j_) 


0  (t) 
n 
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where  E$!t)  is  the  fundamental  matrix  of  (4.63)  o^(t)  is  ^efihed  as  in  (3.45)  and 

<l,i0  t  C([6,“])  is  the  matrix  fulfilling  y^(t)  =  ^(t)(Dg  +  D_)£  ,  where  y^  is  the 

general  C([6,“])-  solution  of  (4.62)  with  f  5  0.  Similarly  a  particular  solution  iji(f) 

of  the  inhomogenous  problem  (4.62)  is  obtained  fulfilling 

-1-6  2max(j0,j_) 

(4.65)  I ) )  ( t)  -  E(Hf )  ( t)  I  <  const,  t  (int)  ap,t> 

where  EHf  is  the  particular  solution  of  (4.3)  which  decays  as  fast  as  possible  according 


to  Theorem  4.3  and  (4.58).  o_  is  defined  as: 

P 


(4.66) 


I (Hf ) <  t) I  -  0(0  ( t) )  . 

P 


Theorem  3.2  remains  valid  if  E$(1)(D  +  d_)  is  substituted  by  i|i°(1)  where  has  been 

continued  to  [I,011]  . 

If  only  an  existence  theorem  for  <J^(t)  is  desired  without  taking  care  of  the  decay 

properties  then  i  =  r  and  8  >  -1  in  (4.61)  can  be  chosen.  In  this  case  the  right  hand 

sides  of  (4.64),  (4.65)  do  not  contain  eventually  occurring  exponentially  decaying  factors 

in  a.  (t)  reap,  a  (t)  . 
h  p 

Moreover  it  is  important  to  consider  problems  where  the  matrix  A(  t)  defined  in 
(4.61)  decays  exponentially: 

(4.67) 


A(t)  =  a(t)tTeq(t*  ,  a  e  Cfa{  [  1 . “» ) ) ,  q(t)  ♦ 


Using  the  same  methods  as  in  the  case  of  algebraic  decay  we  get: 

(4.68)  l*°(t)  -  E*(t)(D0  +  D_)l  <  const.  tY+(ci+1  ^(int)  0  "  oh(t)eq(t) 


v+rn+iii  2max( j  ,  j  )  ,  . 

(4.69)  l(i|i  (f))(t)  -  E  (Hf )  ( t)  I  4  const.  tT  '  tint)  0  (t)eqlC'  . 

P  P 

Only  the  construction  of  the  'particular*  solutions  has  to  be  changed  in  order  to  get  these 

estimates. 

Finally  the  author  conjectures  that  it  is  possible  to  change  r  to  r  defined  in 
Chapter  2  so  that  all  statements  made  should  hold  with  r  instead  of  r  .  Maybe  someone 
else  is  successful  in  improving  the  derived  estimates. 
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5.  Nonlinear  Problems. 


We  consider  problems  of  the  form 


y*  “  t  f  ( t,y) ,  1  <  t  <  •»  ,  a  e  NQ 

y  e  C( [1,-1  ) 
b(y(1 ),  y(“))  ■  0  . 


Prom  (5.2)  and  (5.1)  we  conclude  that 


Now  we  define  for  a  £  R  ,  x,t  e  R  : 


f(»,y(®))  =  0 


and  assume  that 


Sx(a>  -  {y  £  R  |ly  -  al  <  x} 

C  (t,a)  =  { ( t,y )  £  Rn+1|t  >  t,  y  £  S  (a)} 


£'fy  £  Clip(Cx(1'y(*,» 


for  some  x  >  0  sufficiently  large. 

(5.4)  is  a  (nonlinear)  system  of  equations  from  which  y(“)  can  be  calculated  as  a 

"l 

solution  manifold  yC”)  -  y^fp)  ,  yeScR  n^  <  n  ,  if  the  Problems  (5.1),  (5.2), 
(5.3)  admits  a  solution.  The  dimension  of  this  manifold  -n^-  can  be  calculated  a  priori 
if  we  require  for  a  solution  point  y^  that 

(5.8)  ranktf^-.y^))  -  n  -  ^ 

(5.9)  rank(f^(“,y) )  <  n  -  n1  for  y  £  S^y^),  x  >  0  . 

Then  there  is  a  x^  >  0  and  a  ^-dimensional  manifold  y^lp)  which  fulfills  the  equation 

"l 

(5.10)  ff-.y^Ui))  HO  for  u  £  S  c  r 

(5.11)  yjM)  e  S  (y  )  • 

1 

From  (5.8)  we  conclude  that  n1  equals  the  geometrical  multiplicity  of  the  eigenvalue  zero 
of  the  matrix  f  (“,y m)  . 

However  as  practical  examples  point  out,  the  assumption  that  f  («*,y)  does  not 
decrease  its  rank  in  ya  is  too  strong  and  therefore  we  regard  n^  as  a  priori  unknown  but 
obviously  the  solution  of  the  equation  (5.1)  determines  n1  for  a  practical  problem. 
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Now  we  define 


(5.12)  A(t,y)  -  fy(t,y JM)) 
and  require  that  A  fulfills  (3.4),  (3.5)  so  that 

OO 

(5.13)  A(t,u)  =  l  A  (u)t_i  for  t  >  6  . 

i  -  0 

We  transform  Aq  ( *j  >  to  its  Jordan  canonical  form  (y  ) 

(5.14)  AQ(y)  =  E(U)J0(|i)E_1  (y) 
and  introduce  z  as  new  dependent  variable 

(5.15)  E(y)z  =  y  -  y^fy) 
getting 

(5.16)  z'  “  taj(t,u)z  +  tag(z,t,li) 

(5.17)  z(-)  =  0  . 

Here 

(5.18)  J(t,y)  =  E-1()l)A{t,U)E(U> 
and 

(5.19)  g(z,t,U)  ”  E  1  (li)f(t,E(u)z  +  yw(U))  -  J(t,y)z 
hold.  The  perturbation  g  fulfills  the  estimates: 

(5.20)  lg(z,t,M)l  <  C1(M)(lf(t,yo>(U))l  +  »zl2) 

(5.21)  Igfz^t.y)  -  g(z2,t,y)l  <  C2<M>(,Z1I  +  ,z2,),zi  “  z2' 
where  C^(y)  depend  also  on  the  Lipschitz-constant  of  f  on  Cx( 1 ry^lu ) ) . 

We  restrict  y  to  subsets  S  c  s  which  are  defined  as  follows. 

1)  The  projections  onto  the  direct  sums  of  invariant  subspaces  of  JQ(y)  »  which  belong 
to  eigenvalues  with  positive,  zero  and  negative  real  part  are  constant  for  u  t  S. 

Moreover  the  projections  onto  the  invariant  subspaces  of  ( u )  are  constant  for  yes, 
therefore  r+,rQ,r_,r  are  defined  for  JQ(y)  as  in  the  last  chapters  and  are  independent 
of  yes. 

2)  y^fy),  E(y),  E  1 (p )  are  continuous  for  yes. 


Let  $(t,p)  be  the  fundamental  matrix  of  the  (homogenous)  Problem  (5.16). 
3)  The  same  columns  vj(t,p)  of  the  matrix  $(t,P)  fulfill 


~ (a+1 ) r-e ( p)  j. 

I^Mt.p)!  <  CJpJt  (£nt)  1,  e^p)  >  0 

for  all  p  e  S.  Therefore  there  is  a  projection  matrix  independent  of  p  in  S  so 


that 

(5.23) 


l*(t,p)(D0  +  DJI  <  C(p)t 


-(a+llr-E^P) 


(£nt)  for  p  £  S 


holds.  rQ  be  the  number  of  I's  in  the  main  diagonal  of  . 


We  require  f  to  fulfill 

(5.24) 

For  p  £  S,  C  £  C*‘  fixed  we  set 

(5.25) 

where  H  is  as  in  (4.56)  (with  E  =  I).  Regarding  i|i(*,p)  as  an  operator  on  the  Banach 
space : 


-2  (a+1 )  r-£  (p  ) 

Iftt.yjp))!  <  C(P)t  ,  E2(P)  >  0  for  P  £  S 


(♦(z,P))(t)  =  <t>  ( t,p )  (Dq  +  D_)?  +  (Hg(  z,*  >p  ) )  (t) 


A  .  -  { z  |  z  { t )  -  Z(t)t"(“+1,r'e(£nt)  °,  Z  £  C.  ([«,“>)},  Izl  =  IZ*  r  A  .  ), 
t,6  b  £  [6,“>] 


(5.26) 

6  >  1,  0  <  £  =  mini  £1  ,e2>  . 

Every  fixed  point  z  of  \|>(*»P)  establishes  a  solution  for  all  £  £  C  .  At  first  t 
maps  A£  ^  on  A£  j  for  6  sufficiently  large  because  of  (5.20),  (5.23),  (5.24),  (4.55), 
(4.56)  . 

Now  we  take  a  sufficiently  large  sphere  S£(p)  in  A£  ^  with  center 
$(*,p)(D  +  Dj£  and  radius  p  and  prove  the  contraction  property  of  ij)  on  SJ6) 

•  ♦(z^p)  -  Kij,n)lE  -  lH(g(z^  ,  •  ,p )  -  g(  z2  ,•  ,p  ) )  «e  < 

(5.27) 

-e  230 

<  const(p)  P’S  (InS)  Iz^  -  z^l £ 

for  zy'z2  E  Sjp)  because  of  (5.21),  (4.55)  and  (4.56). 

Moreover  if  z  e  SJp)  then 


(5.28) 


l*(z,u)  -  ♦  (•,p)(D  +  D  ) £. * E  “  IHg(z,*,p)l 

2  -£  2  30 
<  (const(p)+p)  •  6  (inS) 


-36- 


Therefore  I (i(r,w)  E  S  (P )  if  5  is  sufficiently  large  and  from  (5.27),  (5.28)  we  conclude 
that  has  a  unique  fixed  point  z  £  S£(p)  A£  ^  for  6  sufficiently  large.  The 

construction  of  H  implies  that  (DQ  +  D  )P(6,p)  ^(O)  *  <5D*^  »W  >  (D^  +  D  ) £ . 

Therefore,  for  fixed  w  e  s  ,  we  have  constructed  a  rQ  +  r  -dimensional  solution 
manifold  in  Re  ^  '  for  ®  sufficiently  large  but  fixed  whenever  ( ( DQ  +  D  K  varies 

rQ  +  r_ 

in  a  compact  set  K  c  C  .  In  order  to  get  more  information  on  the  asymptotic 

behavior  of  the  solution  we  now  treat  the  important  case: 

(5.29)  f ( t , y^ ( w ) )  SO  for  t  >  S  ,  u  £  S 
and 

(5.30)  1 4>  ( t ,  u  )  (  Dq  +  D_)l  ■=  plt^le^^'^t^1*^  (tnt)^,  p  £  C([1,“]) 
where 

(5.31)  q(  t ,(j  )  +  -<•  for  t  +  “  and  u  £  S 
holds. 


We  define: 
(5.32) 


o(t,p> 


eq(t,U)t0(U)(4nt)jo 


and  set 

(5.33)  ($(z,U))(t)  -  ♦(t,u)(D0  +  D_)C  +  (Hg(  z,  •  ,(i )  )  ( t)  ,  (  £  c"  . 

H  is  constructed  according  to  Chapter  4  (with  E  =  I)  so  that  inhomogenities  f  which 
2 

decay  as  o  (t,p)  produce  a  particular  solution  Hf  which  decays  as 
(q+1  )5  2  »>ax(j0,j_) 

t  a  (t,y)(£nt)  .  We  regard  :|i  as  an  operator  on  the  Banach  space 

(5.34)  (A  ,  -  { u ! u  =  o(t,U)U,  U  £  C.(16,»))},  lul  =  101 . .  ,)  . 

The  contraction  mapping  theorem,  employed  as  in  the  case  of  algbraic  decay,  assures  the 
existence  of  a  unique  fixed  point  z  in  j  . 

From  (5.33)  we  conclude 

(5.35)  I  z(  t)  -  ♦<t,u)(D()  +  D_)5»  4  C(u)t(a+1,V(t,lJ)(£nt);' 
where  j  -  max(j0,j_). 

It  is  easy  to  check  that  <H',U)  is  also  a  contraction  in  a  sphere  around 
♦  C.PMDjj  +  D_)£  in  Ao  6  . 


The  uniqueness  of  the  fixed  point  assures  that 


for  every  fixed  point  z  of  '('(•/U)  in  A  ,  .  Therefore 

€  ,  0 

(5.37)  (*(z,|i))(t)  »  $(t,u)(D0  +  D_)C  +  ( (H  -  H)g(z,* ,p) )  (t)  +  (Hg(z,* ,u ) ) ( t ) 

holds.  Because  Hg  and  Hg  are  particular  solutions  for  fixed  z  £  A  ,  we  get 

0,0 

(5.38)  ( (H  -  H)g(z,*,u)(t)  «  $(t,u)(D0  +  d_)y(z),  y  ( z)  e  cn. 

Choosing 

(5.39)  5  =  C  +  Y  ( z) 

assures  that  every  fixed  point  of  t|i(z,u)  is  also  a  fixed  point  of  ijl(z,p)  (with 
C  different  to  £)  and  vice  versa. 

In  general  our  perturbation  approach  does  not  give  us  all  C(  1 1 ,°») )  solutions  of  the 
Problem  (5.1),  (5.2),  (5.3),  it  only  gives  us  all  solutions,  which  decay  at  least  as  fast 
as  t  ^a+1^r  e  where  e  >  0  .  This  is  illustrated  by  the  example: 

(5.40)  y'  =  y2  ,  1  <  t  <  » 

(5.41)  y  £  C([1,»)) 

which  has  the  solution  manifold  y  =  -  -  ^  g  ,  c  >  -1  and  y  i  0  .  Our  approach  gives 

(5.42)  y„  -  0  ,  f^yj  -  0  ,  r  -  i 
and  therefore  the  fixed  point  equation 

(5.43)  y  =  Hy2 

results.  In  A  .  the  only  solution  of  (5.43)  is  y  =  0  ,  which  is  the  only  solution  of 
(5.40)  decaying  faster  than  t-1  . 

If  f  (t.y^tu))  is  not  analytical  at  t  =  •  but  if  it  fulfills  a  relation  of  the 
0 

form  (4.61)  then  i)i  (t,p)(D  +  )  resp.  ij'(  g(  z^  ,•  ,)i ) )  defined  in  Chapter  3.4  have  to  be 

used  instead  of  E(li  )$  ( t,u )  (Dg  +  D  )  resp.  E( U  )Hg(  z, •  ,U ) .  The  results  do  not  change. 

The  following  theorem  follows  immediately: 

Theorem  5.1.  Let  f  fulfill  a  uniform  Lipschitz  condition  in  y  on  [1,"].  Then  the 

problem  (5.1),  (5.2),  (5.3)  where  (5.13)  holds  asymptotically  has  a  solution 
y  i  y(  •  ,  (Dg  +  D_)5,P)  =  E(u)z(*,(DQ  +  D_)5,u)  +  y^iM)  for  every  root  ( ( DQ  +  D_ )  C ,  W  ) 


of  the  equation 


b(E(u)*O,(D0  +  D_)5,U)  +  y— tw) ,y^(w) )  -  0 


_  ^  j.  “RT«^T**4  A 

where  <C,w)  e  C  *  S  and  biR  +  R  .  Here  s(t,(DQ  ■*■  d_)£,u)  denotes  the 

continuated  fixed  points  of  V>(*,W)  with  tec".  On  the  other  hand,  if  the  boundary 

value  problem  (5.1),  (5.2),  (5.3)  has  a  solution  y  ,  so  that  y  -  y(«)  e  A  for  some 

_  £ »  » 

A  n 

£  >  0  ,  then  there  is  a  u  £  R  and  a  (0Q  +  D_)£  ,  C  £  c  so  that  z  :* 

E  1(n)(y  -  y  )  has  the  asymptotic  expansion  (5.35). 


f 


6.  Cage  Studies. 

In  this  chapter  problems  from  fluid  dynamics  reap,  thermo  dynamics  described  by 

boundary  value  problems  on  infinite  intervals  are  analyzed*  These  problems  are  represnted 

by  nonautonomous  system  of  nonlinear  differential  equations*  An  autonomous  problem# 

namely  von  Karman's  swirling  flow  problem  has  been  investigated  by  Markovich  (19R0b).  The 

following  equation  represents  a  model  for  viscous  flow  past  a  solid  at  low  Reynolds  number: 

k  2 

(6.1)  u"  +  — u1  +  ouu'  +  8(u‘ )  “0,  i  4  x  <  "  ,  a  >  0 

(6.2)  u(8)  «  0  ,  6  >  0  ,  u(»)  -  U  >  0  . 

The  parameters  and  variables  are  described  by  Cohen,  Pokers  and  Lagerstrom  (1978).  The 
transformation 

(6-3)  y,  -  n  ,  y2  -  u'  ,  y  -  (y1(y2)T 

takes  (6.1),  (6.2)  into 


(6.4) 

(6  5) 
(6.6) 


y' 


-  y2  -  -  ey2 


f(x,y) 


(a)  (1,0]y(«)  -  0  ,  (b)  [1 ,0]y(»)  «  U 

y  c  C([5,«n  . 


r  O  i 

Ws  get  ym  »  ym(U)  “  l  0  J  and  calculate 


(6.7) 


Vx'y-<u,) 


0 

1 

♦ 1 

0 

0 

0 

-au 

X 

0 

-k 

A0(O) 


A0(U)  has  the  distinct  eigenvalues  0  and  -au  and  we  get 


(6.8) 


v°) 


0  0 

r  i  .  i-i 

,  E(U)  - 

aU 

0  -au 

0  1 

—  — ' 

_  - 

The  transformation  (5.15)  with  u  =  U  results  in  the  system 


0  0  ,  0  - - 

1  aO 

+  — 

.0  -au  X  ,  0  -k 


z  +  g(z,U) 


z<“)  =  0 


The  homogenous  problem  z'  «  (J  (U)  +  —  J  (U))z.  has  a  fundamental  matrix  <)>(x,U)  of  the 

fj  u  xi  n 


<Mx,U)  -  P(x,U) 


10  10 

.  -aUx  .  -k 

0  e  Ox 


(6.12)  P(x,U)  =  I  +  0(x  ') 

holds.  This  follows  from  Theorem  3.1  and  from  the  algorithm  for  the  calculation  of  the 


coefficients.  Moreover 


D0  -  0  .  D.  = 


holds  and  the  fixed  point  equation 


z ( t )  «  P(t,U> 


-oUx  -k_ 
e  x  5 


+  (Hg(z,U) )  ( t) ,  £  e  R 


results  because  we  are  only  interested  in  real  solutions. 
g(z,U)  fulfills 


(6.15) 


lg!z,U)l  <  const(U)lzl 


because  ftx.y^lU))  -  0  holds. 

Chapter  5  assures  the  existence  of  solutions  z(*,5,U)  eh  = 
-1-e  - 

•  {u  |  u  ( x )  -  x  JUixU(x),  U  e  C.  (15,“) )  and  from  (5.35)  we  conclude: 


(6.16) 


Iz(x.C.O)  -  P(x,0) 


-aux  -k_ 
e  x  C 


.  1_2’'  "2aUX  0 

I  <  C(U)x  e  (tnx)2 


Using  Theorem  3.4  this  estimate  can  be  improved,  so  that  its  right  hand  side  is  of  the 


-2k  -2aUx  _  .  .. 

order  x  e  .  Resubstituting  we  get 
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,  ....  „  1  -1.,  -cxtlx  -k,  .  -2aUx  -2k. 

(6.17)  u(x,5,u)  »  u - -  (1  +  0( x  ))e  x  5  +  0(e  x  )  ,  x  ♦  ■  . 

au 

Since  O  e  R+  is  given  (6.5) (a)  has  to  be  used  in  order  to  determine  (  c  K  and  the 
problem  is  well  posed  concerning  the  number  of  conditions  as  x  «  6  and  x  “  “ 

The  second  problem  is  a  similarity  equation  for  a  combined  forced  and  free  convection 
flow  over  a  horizontal  plate  (see  Schneider  (1979)): 


2f».  +  ff»  +  yxg  m  o 
2g‘  +  fg  -  0 


0  <  x  < 


(6.19)  (a)  f(0)  -  f'(0)  -  0  ,  g(0 )  -  1  ,  (b)  f*  («)  -  1  ,  9(“)  =  0 

The  variables  and  the  parameter  k  are  explained  in  Schneider  (1979).  For  simplicity  we 
have  set  the  Prandtl  number  to  1  .  Because  of  (6.19)(b)  we  substitute 


f(x)  «x  +h(x) 


and  get  the  new  problems: 


(a)  2h" '  +  (x  +  h)h"  kxg  «  0 


(b)  2g'  +  (x  ♦  h)g  -  0 


0  <  x  <  “ 


(6.22)  (a)  h( 0 )  -  0  ,  h'(0)  -  -1  ,  g(0)  -  1  ,  (b)  h'(”)  -  0  ,  g(»)  *  0  . 


Substituting 


yi  ”  h  '  y2  “  h<  '  y3  “  h"  *  y4  “  q  '  y  ”  (yl'y2'y3'y4> 


we  get  the  system 


(6.24) 


y‘ 


.  1 


+  5T)y3  -  2  y4 


_  1 


'/2(1  +—  >y4 


xf(x,y)  ,  0  < 


(6.25)(a) 


- 

— 

1 

0 

0 

0 

i 

0 

1 

0 

0 

y(0)  - 

-i 

,  (b) 

0 

0 

0 

1 

1 

_ 

_ 

- 

0  10  0 


0  0  0  1 


We  only  admit  solutions  fulfilling  h(“)  e  R  therefore  we  reouire  that 

(6.26)  yeC(t0,-)> 
holds.  From  (6.23),  (6.24)  we  conclude 

(6.27)  y  -  y  ( h  )  -  (h  ,0,0,0)  ,  h  -  h(-)  e  R  . 

00  00  00  00  oo 

We  calculate: 


(6.28) 


Vx,y-(h-)) 


0  0  0 
0  0  0 


0  -  V2  -  I 


o  -  V2 

_ J  ' 


♦  ! 
X 


0  *r 


0 


A,(hJ 


and 

(6.29) 


E  =  E( h^)  »  diag(1,1,1,-  ^  )  • 


The  substitution  E(hw)z  -  y  -  y^thj  gives  the  system 


(6.30) 


z' 


r 

0 

0 

0 

0 

i _ 

— \ 
0 

0 

0 

_ J 

«  X 

0000 

1 

+  — 

0010 

h 

0 

0 

1 

X 

0  0  “  0 

h 

00  0  -  v2 

1  J  * 

< - 

.  O 

O 

O 

■ 

1 _ 

z  +  xg(z,h  ) 


J0  Vh-)  '  W 

Jq  has  the  eigenvalues  0  with  algebraic  and  geometric  multiplicity  2  and  -  V2  with 
algebraic  multiplicity  2  and  geometric  multiplicity  1  •  Because  f (x.y^lh^) )  «  0  holds 
we  get 

(6.31 )  (a)  lg(z,hjl  4  Cthjlzl2. 

We  have  to  set  up  the  fundamental  matrix  $ ( x , h^ )  of  the  system 

(6.31 )  (b)  X(J0  +  ^  J1(h-,,Zh  '  2h=  (V2h'2h'Zh’- 


Because  of  the  simple  structure  of  we  do  not  have  to  apply  the  algorithm  of 

Chapter  4  ,  we  can  proceed  in  the  following  way. 

The  last  equation  of  (6.31)  is 


(6.32) 


4'  .  1.  -1,4 

rh  -  x(-  V2-  -)*h  • 


It  can  be  integrated  at  once  giving 


(6.33) 


-  V4xZ  -  x 
e  C 


C  e  R 


Setting  c  -  1  we  have  to  find  a  particular  solution  of 

h 


(6.34) 

We  take 

(6.35) 
Integrating 

(6.36) 


-  *<- V  r*  ) 


3  ,  4 

,  ,z.  +  xz.  . 

2  x  h  h 


1/  2  -1 

-  74  x - -  —  2 

*  2  x  x 

6  —  ( 1 - - 

2  2 
x 


1 
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we  find 


(6.37) 


1.2  “ 

2  "  '/4x  -  —  x 

z  «  e  x  0(1) 


Analogously  we  integrate 

(6.38) 

and  get 


1'  2 
Zh  “  *h 


(6.39) 


1/  2  h" 
i  -  4*  ~T 


0(1) 


The  fourth  column  of  can  ^  chogen  afl  <VVV*4)T. 

In  order  to  get  the  third  column  we  set  C  «  0  in  (6.33)  and  proceed  as  we  did. 
Finally  we  get  the  following  fundamental  matrix 


r 


(6.40)  $(x,h  ) 


1  x 


0  1  e 


0  0  e 


0  0 


2  h 
x  " 

4  "  2 


2  h 


0(x'2)  e 


4  2 


2  h 
x  «° 

4  ”  2 


2  h 
x  “ 


0(1) 


4  2  0(x-1)  e  4  2  x  01 1 ) 


2 

X  00 

”  4  ”  2  x  2,,  1  , 

®  X  d - 2  ) 

X 

2  h 

x  00 

-  4-  -  r  x 


which  we  can  write  as  in  Theorem  4. 1 


Q(x,h  ) 


(6.41) 

where 


♦  ( *  *  h  )  »  P(x,h  )x  e 
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1 


1 


0 


0 


(6.42) 


P(x,h_) 


0  0 
0  0 
0  0 


+  -0,1, 


(6.43) 


D  »  diag(1,x,1,x  ) 


(6.44) 
holds . 


2  h  2  h 

x  “  x  • 

'  - - T—  X  -  —  -  —  X 

Qtx.h^)  -  diag(l,1,e  ,  e  ) 


Slmllarily  to  the  first  example  we  conclude  that  there  are  solutions  *(x,C..C_,h  ) 

i  2  00 

in  A  which  is  now  the  space  of  all  functions  in  C( [{,“))  which  decay  at  least  as 

e ,  o 

x  *  £ inx  ,  e  >  0  because  r  =  2  and  a  =  1  holds.  Moreover 


(6.45) 


xCx.C^Cj  *»_)  -  Pfx.h^J 


0 

0 

x2  h- 

4  "  2  X  , 

xe  5, 


2  h 
x  » 

2  ‘  4  "  2  * 

x  e  5 


2  -I 


-  -nr  -  h  x 
,  4  2  ”, 

+  0(x  e  ) 


with  h^.S^Cj  R  results  because  the  exponential  factor  e 

the  fundamental  matrix. 

So  we  get  the  asymptotic  expansions 

2  h  2 

x  “  x  , 

- - — —  x  -  — h^x 

(6.46)  f(x)  -  x  +  h  +  0(1 )e  4  2  +  0( xe  ”  ) 


~  '  h»x 


does  not  appear  in 


2  h  2 

x  ■  x  . 

2  _  4  2  X  4_2~ 

(6.47)  g(x)  -  -  £  e  *  C2  +  0(x  e  ) 

where  the  0(1)  in  (6.46)  depends  linearly  on  £  and 
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The  constants  <h„/C , ,?2 )  e  R3  have  to  be  determined  from  the  three  initial 
conditions  (6.2.’i)(a). 

The  third  problem  to  be  analyze d  is  the  well  known  Falkner-Skan  equation  laee  for 
example  Lentini  and  Pereyra  ( 1977 ) ) : 

(6.48)  f"'  ff"  +  (1  -  f'2)  -  0 

(6.49)  f'(«)  -  1. 

We  do  not  pose  any  initial  condition  because  we  look  for  a  solution  manifold. 

Because  of  (6.48)  we  substitute 

(6.50)  f(x)  -  X  +  g(x)  ,  y,  -  q  ,  y2  -  g‘  ,  y3  -  g" 


T 

(6.51)  y  -  (y1#  y2,  y3) 
and  get  the  system 


(6.52)  y'  -  x 


r 

Moreover  we  require 

(6.53) 


y  c  C( [6 ,»]  ) 


-  x  f(x,y) ,  x  *  6  , 


so  that 


(6.54) 

holds,  we  calculate 


y(")  -  <g„,0(0)T  , 


9.  *"  g(") 


(6.55) 


Vx'Y-(q-)) 


The  linearized  system  is  therefore 


(6.56)  z'  «  x(J  +  —  J,(g„))t. 

fl  U  X  1  n 

Because  Jg  has  the  eigenvalue  0  with  algebraic  multiplicity  2  we  apply  the  theory 
developed  in  Chapter  4. 
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At  first  we  split  up  the  system  (6.56)  by  the  transformation 


(6.57)  *h  ■  P  ( x , g^ ) u  u  -  (u1(u2,u3)T,  Ptx.g^)  ~  I  +  £  P^tg^lx  1 

to  get  subsystems  whose  leading  matrices  have  the  only  eigenvalue  0  reap.  -1.  (6.57) 

gives  a  system  of  the  form; 

(6.58)  u*  -  x  B(x,g<o)u  ,  Bfx.g^)  ~  JQ  +  £  B^fg^Jx  1  . 


i”1 


From  Wasow  (1965)  we  conclude  that 


0 

0 

Pil" 

'bi, 

bi2 

0 

(6.59) 

Pi  ' 

0 

0 

Pi2 

'  Bi  " 

bi3 

bi4 

0 

.  Pi3 

Pi3 

0 

0 

0 

bi5  . 

holds  and 

that  the  recursion 

i-1 

(6.60) 

Jo  pi  - 

Vo  " 

l  <p 
-  8 
s»0 

Bi-s 

-  Ji-s  Ps>  - 

a  - 

2,Pi-2  ' 

i  >  0 

with  the 

last  term  absent  for 

1  <  2 

and 

Jy  »  0  for 

k  > 

1  holds. 

From  the  investigation  of  the  perturbed  system  we  know  that  only  the  coefficients  B0, 
b^.Bj  ,  B-j  and  B^  influence  the  asymptotic  behavior  of  the  fundamental  matrix  because 
(a  +  1)r  «  4  for  our  example.  (6.59)  and  (6.60)  give 


- 

- 

0  1  0 

0  0  0 

0  0  0 

'  °2  " 

0  2  0 

o 

O 

1 

9° 

0  0-2 

j 

- 

(6.61 ) 


"  0 

0 

0 

-  0 

0 

0  - 

B3  - 

0 

_  0 

-Z'J. 

0 

0 

-29.  _ 

B4  “ 

0 

0 

-2+2<*i 

0 

0 

b45  - 

We  do  not  have  to  know  b^j  explicitely  because  it  does  not  influence  the  behavior  of  the 
solution  of 


(6.62) 


+  •••  >u3  . 
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From  Chapter  3  we  get 


1  '  g«X  -2 


3  r3' 

Moreover  we  get: 


'  3' 


■1  ,  .  -2 
^3 1 '  +  P32 


(6.64) 


:m:  ;]•?[: :]•?[:  jj 


n  .  +  * •• J v 
0  -2+2gf  J  U2 


•)(„ 


because  the  leading  terra  which  comes  from  Bg  vanishes.  Therefore  the  coefficient  of 
does  not  influence  the  behavior  of  u^Uj  *  is  sufficient  to  9olve 


(6.65) 


0*M)  -°J 


We  get  directly  by  integration 


(6.66) 


2g„ 


U2  “  X 


exp(— )  and  =  x  0(1) 


Finally  we  conclude 


(6.67)  f  *  x  +  g  +  p(x,g  ]« 


g  x  -x  -2g  x 

OO  /  3od 


5  +  o(e 


),  Pfx.gJ  -  1  +  l  Pi<9jX 
i-1 


-i 


with  5  £  R  because  we  look  for  real  solutions  and  because  there  is  no  column  of  the 


fundamental  matrix  of  (6.56)  which  has  the  factor  e 


-x‘  -2g  x 
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